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First-Principles Study of Electronic Structure in a-(BEDT-TTF)2l3 
at Ambient Pressure and with Uniaxial Strain. 
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Within the framework of the density functional theory, we calculate the electronic structure of 
a-(BEDT-TTF)2l3 at 8 K and room temperature at ambient pressure and with uniaxial strain 
along the a- and &-axes. We confirm the existence of anisotropic Dirac cone dispersion near the 
chemical potential. We also extract the orthogonal tight-binding parameters to analyze physical 
properties. An investigation of the electronic structure near the chemical potential clarifies that 
effects of uniaxial strain along the a-axis is different from that along the 6-axis. The carrier 
densities show dependence at low temperatures, which may explain the experimental findings 
not only qualitatively but also quantitatively. 
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1. Introduction 

A quasi-two-dimensional organic charge transfer con- 
ductor, Q;-(BEDT-TTF)2l3, has attracted much atten- 
tion due to its exotic physical properties. The tempera- 
ture dependence of the resistivity shows a metal-insulator 
I transition at 135 K, which is induced by charge-ordering 

■ at low temperatures. When hydrostatic pressure or 
' uniaxial strain suppresses the charge-ordering, the tem- 
' perature dependence of the resistivity shows a metallic 
] behavior even at low temperatures. However, the trans- 

■ port property in this metallic state is very unique. It is 
[ reported that the carrier density decreases and the mo- 

■ bility increases drastically with decreasing temperature; 
I this results in temperature- independent resistivity.^ Ka- 
' jita et al. named this state the narrow-gap semiconduct- 
I ing state. 

Applying uniaxial strain is now one of the important 
] techniques of changing the physical properties of organic 

■ charge transfer salts. ^ Tajima et al. reported that the 
transport property depends on the direction of uniaxial 
strain. They discovered superconductivity with 1-5 kbar 

I of uniaxial strain along the a-axis.^ However, no super- 

■ conductivity appears with uniaxial strain along the b- 
[ axis. In addition to this, the metal-insulator transition 

■ temperature decreases when uniaxial strain is applied 
along the a-axis, while the metal-insulator transition 
temperature only nominally changes with uniaxial strain 
along the 6-axis. Thus the obtained phase diagrams are 
different depending on the directions of uniaxial strain. 
Kondo and Kagoshima measured the change of the crys- 
tal structures with uniaxial strain at 8 K and room tem- 
perature (RT) to discuss the change of electronic struc- 
tures in the semi-empirical estimations of transfer and 
overlap integrals.^' 

Based on these electronic structures in the semi- 
empirical approximation, Kobayashi et al. discussed the 
charge-ordering, narrow-gap semiconducting and super- 
conducting states.^' ^ Through a careful investigation of 
the electronic structure, Katayama et al. also found a 
unique electronic structure near the Fermi level, called 



two-dimensional anisotropic Dirac cone dispersion, E — 
i\vF ■ k±\, where vp is the Fermi velocity and k± is 
the wave number along the a*6*-plane measured from 
the crossing point. This unique feature may be the ori- 
gin of the narrow-gap semiconducting state, because it 
is qualitatively consistent with the low carrier density 
and compensating high mobility. It also leads us to ex- 
pect large diamagnetism analogous to those observed in 
graphite and bismuth. Although this discovery is very 
important in clarifying the interesting physical proper- 
ties of the title compound, there is a problem that their 
studies rely on the simple semi-empirical approximation. 
Therefore, it is important to confirm that the simplicity 
of their model is not the origin of the interesting elec- 
tronic structure. 

BEDT-TTF charge transfer salts consist of the BEDT- 
TTF conducting and anion insulating layers. The highest 
occupied molecular orbitals (HOMO) of the BEDT-TTF 
molecule usually form the electronic structure near the 
Fermi level. The semi-empirical approximation, which 
evaluates the transfer and overlap integrals between the 
HOMOs of the BEDT-TTF molecules, has achieved great 
success in the prediction of the electronic structure of 
the BEDT-TTF charge transfer organic salts at ambi- 
ent pressure. One of the reasons behind this success is 
that the BEDT-TTF charge transfer organic salts usu- 
ally have large round Fermi surfaces at half- filling, which 
conceals small errors of the estimated transfer and over- 
lap integrals in the extended Hiickel approximation. 

In contrast to the above, the electronic structure of 
the title compound is very delicate; it is semimetallic 
with a very small band overlap (at least in the semi- 
empirical approximation). Though parameters in the 
semi-empirical approximation are tuned up to recon- 
struct the Fermi surface of some of BEDT-TTF com- 
pounds, a small difference in the electronic structure can 
change the Fermi surfaces drastically in the title com- 
pound because of the semimetallicity. Furthermore, the 
electronic structure at high temperatures may be differ- 
ent from that at low temperatures because the thermal 
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contraction reported from the X-ray measurements is ap- 
preciable. Therefore, an accurate first-principles calcula- 
tion is demanded for elucidating the reliable electronic 
structure. In this work, we calculate electronic structures 
for the unit cells at 8 K and at RT, which are determined 
by the X-ray measurements, to clarify the change of their 
electronic structures within the framework of the den- 
sity functional theory and the generalized gradient ap- 
proximation (GGA). In addition, we extract the transfer 
integrals of the tight-binding model to evaluate carrier 
densities and other physical properties. 

We have demonstrated that the GGA has sufficient ac- 
cm-acy to determine the structures of /3'-(BEDT-TTF)2X 
(X=ICl2 and AUCI2) at ambient pressure and various 
appHed pressure. ^-^'^^ Encouraged by these successes, we 
use the GGA also in this study to relax the internal coor- 
dinates and to calculate the electronic structures of the 
title compound at ambient pressure, with uniaxial strain 
of 2 kbar along the a-axis and 3 kbar along the 6-axis. 
We employ the ultrasoft pseudopotential technique with 
planewave basis sets, and the cutoff energies of 25 Ry for 
the wavefunction and 256 Ry for the charge density. For 
A;-point sampling, we use a 4x4x2 mesh during the re- 
laxation of the internal coordinates and an 8x8x4 mesh 
to calculate electronic structures in the relaxed atomic 
positions. In this calculation, we use lattice parameters 
from the X-ray measurements obtained at 8 K and RT, 
because we want to study the thermal effects on the elec- 
tronic structure through the difference in the lattice pa- 
rameters, which is very large in this material. For ex- 
ample, the lattice parameters a and h change by 0.5% 
with 2 kbar strain in the a-axis and by 1% with 3 kbar 
strain in the 6-axis. However, a and b change about 2% 
with decreasing temperature (from RT to 8 K). The pur- 
pose of this study is to evaluate and discuss changes in 
electronic structures depending on temperature and uni- 
axial strain. 

2. Results and Discussion 

2. 1 Results 

2.1.1 Overall electronic structure 

The relaxed internal atomic positions are only slightly 
different from those determined by the X-ray measure- 
ment. The largest difference arises for atoms near the 
terminal hydrogens. In the case of ambient pressure at 
RT and 8 K, the hydrogens move by 0.11 A at most. 
This may be related to the well-known fact that the X- 
ray measurements underestimate C-H bond lengths. The 
other atoms move much less than the hydrogens. Carbon, 
sulfur and iodine atoms move by 0.04, 0.03 and 0.03 A 
at most, respectively. 

We show the electronic structure from -3.5 eV to 3.0 eV 
in Fig. 1. We take the Fermi level to he ai E = 0. 
The HOMO-2, HOMO-1, HOMO and LUMO bands in 
the crystal, which are at -0.5 eV to 0.2eV, are com- 
posed of the HOMOs of the BEDT-TTF molecules. The 
LUMO-M and LUM0-F2 around +leV are composed of 
I;^. There are 14 (7x 2) molecular levels made of two 
below the Fermi level in this energy range. The energy 
bands where the contributions of I3 are large are shown 
by arrows left of the F point, (0,0,0), in Fig. 1. There are 
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Fig. 1. Electronic structure at ambient pressure and 8 K. Arrows 

show the states that have large contributions from 2I3 at the 
F-point, (0,0,0). We take the Fermi level to be at _E = 0. 



more than 14 arrows because mixing occurs between the 
BEDT-TTF and I3 molecules. The bands coming from 
I3" have large dispersion along the a*-axis because I;^ 
molecules form a one-dimensional chain structure along 
the a-axis. 

Let us examine the electronic structure near the Fermi 
level. Figures 2 (a) and 2 (b) show the electronic struc- 
ture at 8 K and RT. The left side shows the electronic 
structure on the a*fe*-plane and the right side shows that 
along the c*-axis. The dispersion along the c*-axis is 
much smaller than those along the other directions, re- 
flecting the two-dimensional nature of the compound. 

Here we compare the bandwidths at ambient pressure 
and with strain. In this paper we define the total band- 
width as the highest energy of the LUMO band minus 
the lowest energy of the HOMO-2 band, because they are 
made of the HOMOs of the BEDT-TTF molecules. The 
total bandwidth is slightly larger with strain of 3 kbar 
along the 6-axis than in the other cases; however, the dif- 
ferences between the a- and 6-axis strain are slight in the 
HOMO and LUMO bands at this energy scale. Figure 2 
(c) shows the electronic structures at ambient pressure 
for 8 K and RT. The changes according to the temper- 
ature are much larger than those according to uniaxial 
strain at the same temperature. Details are summarized 
in Table I. 

Strangely, the bandwidths with strain of 2 kbar along 
the a-axis are almost the same or slightly smaller than 
those at ambient pressure at both 8 K and RT. The band- 
widths with uniaxial strain of 3 kbar along the fe-axis are 
larger than those at ambient pressure. The bandwidths 
of HOMO and LUMO are about 0.1 and 0.13~0.15 eV, 
respectively, at both 8 K and RT. 
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Table I. Bandwidths (eV) calculated in 8x8x4 mesh at ambient 
pressure (ap), with uniaxial strain of 2 kbar along the a-axis (a2) 
and 3 kbar along the fe-axis (b3) at RT (a) and 8 K (b). H-L rep- 
resents the highest energy of the LUMO band minus the lowest 
energy of the HOMO band. Total means the highest energy of 
the LUMO band minus the lowest energy of the HOMO-2 band. 
The state with * is metallic in 8x8x4 mesh. 



T =RT 




ap 


a2* 


b3 


HOMO 


0.098 
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LUMO 


0.136 


0.134 


0.141 


H-L 


0.233 


0.228 


0.213 


total 


o.r)7(i 


0.r)7() 


0.()02 


T =8 K 




ap 


a2 


b3 


HOMO 


0.093 


0.099 


0.099 


LUMO 


0.141 


0.137 


0.151 


H-L 


0.243 


0.243 


0.253 


total 


0.621 


0.616 


0.640 



2.1.2 Anisotropic Dirac cone dispersion 

Next wc concentrate on the 2D electronic structure to 
investigate the existence of the anisotropic Dirac cone 
dispersion discovered by Kobayashi et ah through the 
semi-empirical electronic structure calculation. We dis- 
cuss the dispersion along the interlayer direction later in 
§• 2.3. 

The dispersions of the HOMO and the LUMO from 
{—IT, 0) to (0, tt) appear to cross. We examine this re- 
gion in detail at kz = 0, as shown in Fig. 3. Figure 3 (a) 
shows the electronic structure near (— 7r/2,7r/2) at am- 
bient pressure and RT. It clearly shows the existence of 
the anisotropic Dirac cone dispersion, since the finer we 
make the mesh of the fc-point, the smaller the energy gap 
becomes. The value of the direct gap is 0.54 meV in the 
finest mesh, which corresponds to the 324 x 324 mesh in 
the a*- and 6*-axis. These results strongly support that 
this system has anisotropic Dirac cone dispersion (ac- 
cidental degeneracy) in the a*6*-plane.^^ Note that the 
symmetry of this system is PI and that this accidental 
degeneracy occurs not at a high-symmetry point, but at 
a general point {kx, ky) = (— 0.4947r, 0.5747r). It is known 
that such kinds of degeneracy at the Fermi level take 
place, for example, in graphite sheets and carbon nan- 
otubes at a high-symmetry point. However, the present 
system is unique because the accidental degeneracy oc- 
curs at a general k-point as well as at the Fermi level. 

In Figs. 3 (b) and 3 (c), we project the electronic struc- 
ture in the area from (k^, ky) = (—'K,0) to (0, tt) onto the 
plane perpendicular to the a*&*-plane. Figure 3 (b) shows 
the electronic structures at RT, while Fig. 3 (c) shows 
those at 8 K. Figure 3 (d) schematically represents the 
projection. These figures clearly show the existence of 
the anisotropic Dirac cone dispersion in all cases. At am- 
bient pressure and RT, the energy at (— tt, tt) {E{—tt, tt)) 
is almost the same as that of the crossing point of the 
anisotropic Dirac cone dispersion (i?(cone)). With strain 
of 2 kbar along the a-axis at RT, E{—tt, tt) is lower 
than i?(cone). On the other hand, with strain of 3 kbar 
along the 6-axis at RT, tt, tt) is higher than i?(cone). 
E{—tt,tt) is much higher than E{cone) at 8 K. iJ(cone) 
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Table H. Energy differences (meV) between E(—it, it) and 
i?(cone) at room temperature (RT) and 8 K. The van Hove sin- 
gularity at (— TT, 7r) is the nearest the Fermi energy. i<;(cone) is 
the energy of the crossing point at the Dirac cone dispersion. 
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11.1 


8.6 


6.1 



becomes the Fermi level when it is lower than E{—it, tt) 
if one considers only the plane at kz = 0. Details are 
summarized in Table II. 

There may be errors in the estimated energy. How- 
ever, it is safe to conclude that there is a tendency that 
E{—TT, TT) — E{cone) at 8 K is higher than that at RT, and 
that the Dirac cone dispersion determines the transport 
properties at low temperatures. 

2.2 Discussion 

2.2.1 Bandwidths 

The calculated results show that the total band width 
with uniaxial strain of 3 kbar along the 6-axis is much 
larger than those at ambient pressure and with uniaxial 
strain of 2 kbar along the a-axis. Taking into account 
the fact that the value of strain along the 6-axis is 50% 
larger than that along the a-axis, the change of the total 
bandwidth with strain along the 6-axis is much larger 
(more than 150%) than that with strain along the a-axis. 
Therefore, we can expect that the bandwidth with strain 
along the 6-axis is larger than that with strain along the 
a-axis if the values of strain are the same. The charge- 
ordered phase is usually discussed in terms of transfer 
integrals and long-range Coulomb interaction. Roughly 
speaking, the case where the bandwidth is comparable to 
the Coulomb interaction determines the phase boundary. 
Thus, it is expected that uniaxial strain along the 6-axis 
destroys the CO phase slightly more easily than that with 
strain along the a-axis. This agrees with the experimental 
findings. 

2.2.2 Analysis using tight-binding model for the ab- 
plane 

Next, we extract the parameters of the orthogonal 
tight-binding model to reproduce the first-principles 
band structure in the a6-plane. We have found that the 
transfer integrals of the next nearest neighbor as well as 
those of the nearest neighbor are necessary to reproduce 
the electronic structure in detail. In this paper, we use 
not only the nearest neighbor parameters, al, a2, a3, bl, 
b2, b3 and b4, which have been introduced by Mori et 
al., but also al', a3',and a4', which are the next near- 
est neighbor transfer integrals along the o-axis shown 
m Fig. 4.'^ The parameters arc determined to reproduce 
not only the overall electronic structures, but also the de- 
tailed structure near the Dirac cone dispersion and near 
{—TT, tt). We note that including al', a3' and a4' is not the 
only way to reproduce the band structure. Other sets of 
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Table III. Fitted transfer integrals (eV) at ambient pressure (ap) 
and with uniaxial strain of 2 kbar along the a-axis (a2) and 3 kbar 
along the 6-axis (b3) at room temperature (RT) and 8 K. al-a4, 
bl-b4 and al'-aS' are shown in Fig. 4. al-a4 and bl-b4 under the 
tb column are the transfer integrals used in the semi-empirical 
approximation calculated as overlap integrals times —10 eV and 
are shown for comparison.^ 



T=RT 




ap 


a2 


b3 


tb 


al 


-0.0101 


-0.0173 


-0.0195 


-0.03 


a2 


-0.0476 


-0.0358 


-0.0393 


-0.049 


a3 


0.0093 


0.0029 


0.0287 


0.0018 


bl 


0.1081 


0.1075 


0.1226 


0.123 


b2 


0.1109 


0.1154 


0.1176 


0.142 


b3 


0.0551 


0.0559 


0.0537 


0.062 


b4 


0.0151 


0.0085 


0.0156 


0.023 


al' 


0.0088 


0.0085 


0.0220 




a3' 


0.0019 


0.0072 


-0.0151 




a4' 


0.0009 


-0.0012 


-0.0019 





T=8 K 




ap 
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b3 


al 


-0.0267 


-0.0245 


-0.0267 


a2 


-0.0511 


-0.0461 


-0.0405 


a3 


0.0323 


0.0373 


0.0397 


bl 


0.1241 


0.1276 


0.1336 


b2 


0.1296 


0.1220 


0.1251 


b3 


0.0513 


0.0507 


0.0515 


b4 


0.0152 


0.0133 


0.0190 


al' 


0.0119 


0.0247 


0.0197 


a3' 


0.0046 


-0.0105 


-0.0118 


a4' 


0.0060 


-0.0019 


0.0061 



transfer integrals can also reproduce the electronic struc- 
tures. The Dirac cone dispersion itself can be reproduced 
without al', a3' or a4'. 

We show fitted parameters in Table III. We are not 
able to reproduce relative energies between E{—tt, 0) and 
S(cone) at RT accurately. However, the energy differ- 
ences between the first-principles restilts and the tight- 
binding ones are a few of meV. Such a small difference is 
irrelevant because the thermal broadening of RT is much 
larger. The agreement is much better at 8 K. 

Comparing the transfer integrals at ambient pressure 
and RT with those of the semi-empirical estimation, 
the transfer integrals of bl-b4 are found to agree quite 
well. On the other hand, it would be better to say that 
the transfer integrals along the a-axis cannot give good 
agreement with those in the semi-empirical estimation 
because the next nearest transfer integrals, al',a3' and 
a4', are necessary. The dispersion, as well as the trans- 
fer integrals, along the a*-axis is much smaller than that 
along the &*-axis as shown in Fig. 2. In general, this may 
cause a problem in semi-empirical approximation when 
evahiating smaller values. Therefore, it is not surpris- 
ing that the semi-empirical approximation fails to yield 
precise estimates of transfer integrals along the a-axis. 
These changes of transfer integrals may possibly affect 
the stability of the charge-ordering patterns.^' 15-17 

Here wc compare the first-principles band structure 
with the electronic structure in the somi-cmpirical pa- 
rameters shown in Fig. 57 The overall shapes agree quite 
well. However, we find notable differences in the disper- 
sion near the Fermi energy {E = 0). For example, let 
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US consider the line from (— tt, tt) to (0, tt). The bottom 
of the LUMO band lies at (0,7r) in the semi-empirical 
approximation, but it is at (—■ tt, tt) in the first-principles 
result. In the HOMO band, the top is located at (— tt, 0) 
in the semi-empirical approximation, but 7r/2, 7r/2) 
is higher than tt, 0) in the first-principles result. 
These changes make the Fermi surfaces different from 
each other. This difference in the {—it, tt) — (0, tt) line in 
the semi-empirical approximation is related to the fact 
that the fitted tight-binding parameters along the a-axis 
are quite different from the Hiickel ones. 

The existence of the Dirac cone dispersion gives birth 
to interesting physical properties. We show the density 
of states at each strain and temperature in Fig. 6, using 
the fitted parameters and assuming the dispersion along 
the c*-axis to be irrelevant. E — corresponds to the 
chemical potential in the figure (this includes the effect 
of temperature broadening) . 

First, let us examine the density of states at RT 
(Fig. 6 (a)). There exists a van Hove singularity quite 
near ^^(cone), which makes the density of state at E' = 
finite due to the thermal broadening. Note that the ther- 
mal energy of RT corresponds to about 0.03 eV; therefore 
electronic structures of two peaks above and below E = 
contribute to the physical properties at RT. The DOS 
with strain along the 6-axis looks different from the oth- 
ers. The position of the van Hove singularity above the 
chemical potcintial with uniaxial strain of 3 kbar along 
the 6-axis is higher than those at ambient pressure and 
with uniaxial strain of 2 kbar along the a-axis, as shown 
in Table. II. 

At 8 K (Fig. 6 (b)), the linear behavior of the density 
of states near E = is clear. The nearest peak above 
the chemical potential {E ^ +0.01 eV) corresponds to 
the van Hove singularity at (— tt, tt) and the peak below 
the chemical potential {E ~ —0.02 eV) corresponds to 
another van Hove singularity at (— 7r,0). At this energy 
scale, we can see notable differences among the temper- 
atures and uniaxial strains. While the density of states 
near E = with miiaxial strain of 2 kbar along the a- 
axis is almost the same as that at ambient pressure, the 
density of states with uniaxial strain of 3 kbar along the 
6-axis shows much larger values. 

The property with strain along the 6-axis looks differ- 
ent from the others at both temperatures. The positions 
of van Hove singularities above the chemical potential 
{E{—TT,n) — E{cone) in the first-principles result) are 
summarized in Table II. The peak positions in the den- 
sity of states between RT and 8 K are similar to Ishibashi 
et al.'s result. The difference is mainly due to the fact 
that we take into account of effects of the Dirac cone 
dispersion explicitly. 

Next we calculate effective carrier density as a func- 
tion of temperature. First, we calculate the density of 
states, D{E,T), at each temperature using linearly in- 
terpolated transfer integrals at each strain; then we sim- 
ply estimate the effective carrier densities as n(T) = 
X5j, D(E, T)dE/V{T), where V{T) is the cell volume at 
temperature T. We show the result in Fig. 6 (c). At low 
temperatures, the carrier density is proportional to the 
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square of the temperature, reflecting the hnear density 
of states near the chemical potential shown in Fig. 6 (b). 
At about 100 K, there is a hump. This is due to the peak 
in the density of states above the chemical potential at 
(— TT, tt). This hump would add an electron contribution 
to the Hall effect. The carrier density in our calculation is 
4.8 X 10^^ cm~^ at ambient pressure and 3.6 x 10^^ cm""^ 
at uniaxial strain of 2 kbar along the a-axis at 10 K. The 
experimental carrier density appears to be proportional 
to the square of the temperature with uniaxial strain of 
10 kbar along the a-axis, and its value is about 1 x 10^'' 
cm~^ at 10 K. One cannot compare these results directly, 
because the value of strain is greatly different. However, 
their values agree quite well. 

There is a discrepancy, however, in the sign of the Hall 
coefficient, Rh- Experimentally, the sign of Rh is very 
sensitive. It depends on the value of hydrostatic pressure 
or uniaxial strain. For example, with uniaxial strain of 
10 kbar along the a-axis, R^ is negative at low temper- 
atures, while it is positive at high temperatures. On the 
other hand, with high hydrostatic pressure (20 kbar), it 
is positive over the whole temperature range. ^'^ 

If the particle-hole symmetry is satisfied, Rh = 0. 
Therefore, asymmetry is necessary to give nonzero Rh- 
In the calculated results, the nearest peak from the chem- 
ical potential is located on the positive energy side at 
(— 7r,7r) at 8 K. For example, at ambient pressure or with 
uniaxial strain along the a-axis, E[—Tr, tt) — E(cone) = 
8 ~ 10 meV at 8 K. The second nearest peak from the 
chemical potential is on the negative energy side at (— tt, 
0). However, i?(cone) — £^(—77,0) is twice as large as 
E{—n,Tr) — E{cone) at 8 K. This makes the chemical 
potential slightly lower than £^(cone), which results in 
a positive value of Rh- In contrast, Rh becomes nega- 
tive if E{—TT, 0) (< ii^(cone)) is closer to ^^(cone) than 
£(-77, tt) (> i?(cone)) is. 

Provided that the linear density of states of the Dirac 
cone dispersion governs the temperature dependence of 
the carrier density with uniaxial strain of 10 kbar along 
the a-axis, a possible origin of the change of the sign of 
Rh is that both the peak at (— 7r,7r) on the positive en- 
ergy side and the peak at (— tt, 0) on the negative energy 
side would be easily shifted due to strain. Considering 
that 10 kbar is much larger than 2 kbar, these differences 
could arise. We do not discuss the case of the hydrostatic 
pressure of 20 kbar, because the value is too different 
from strain of 2 or 3 kbar calculated in this study. We 
note that tt, 0) (< i?(cone)) is much closer to the 
^^(cone) than 7r,7r) (> E{cone)) without the struc- 
tural relaxation. 

The experiment also shows that Rh > at higher 
temperatures. However, it is difficult to conclude which 
contribution is dominant at RT because the dispersion is 
very complex. The experimental crystal structures with 
high hydrostatic pressure and with high uniaxial strain 
are unknown. Further theoretical and experimental stud- 
ies on this problem are necessary. 

In this paper, we do not discuss the origin, or the con- 
dition of the existence, of the anisotropic Dirac cone dis- 
persion at the Fermi level and at a general k-point. Here, 
we mention that these first-principles calculations and 
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their tight-binding fitting show that the situation is ex- 
pected to be realized in a rather wide range of parameters 
(temperature and strain). The condition of the existence 
in simple tight-binding models can be found in ref. 10. 

2-3 Electronic structure along the c-axis 

Next, we consider the dispersion along the c*-axis, 
which is thought to be very small and has so far been ne- 
glected in the semi-empirical parameterization and tight- 
binding models. However, there is a possibility that the 
transfer integrals along the c-axis are relevant at low tem- 
peratures because of the semimetallic nature. 

First let us examine the overall nature. Table IV shows 
the energy differences, E{kz = 0) — E{kz = tt), along the 
c*-axis. (fcj;, fcj,)=(— 7r/2,7r/2) is near the crossing point 
of the Dirac cone dispersion. The band dispersion along 
the c*-axis is 4 meV at the greatest. The band disper- 
sion of the HOMO band at (— 7r,0) is very small. These 
reflect the fact that the nearest molecule perpendicular 
to the a6-plane is along the (1/2,1/2,1) direction. These 
values are clearly much smaller than those of the band 
dispersion in the a*/j*-plane. 

It has been thought that this material is a two- 
dimensional conductor because the resistivity along the 
c-axis is 10^ times larger than those in the a6-plane 
(at ambient pressure). Whether the transfer integrals 
along the c*-axis (t^) are relevant or not is an interest- 
ing problem. If carriers are confined in the conducting 
layer, or conductance along the c-axis does not occur via 
coherent hopping, then it will be safe to neglect t^ in 
the analysis of transport properties within the a6-plane. 
However, may be relevant at temperatures lower than 
the energy scale of tz ~ 10 K; then we must take into 
account of the effect of t^- If the coherent hopping occurs 
between BEDT-TTF molecules directly along the c-axis, 
more calculations are necessary. 

We briefiy look at the dispersion near the Dirac cone 
dispersion. We show the dispersion along the c*-axis at 
some fc-point near the crossing point of the Dirac cone 
dispersion in Fig. 7. It is not easy to predict the exact 
fc-point where the gap is zero because this is a numerical 
study. The figm-e shows that a gapless point exists at 
kz ~ O.Stt when the energy gap at fc^ = is not zero. We 
expect that the gapless kz-po'mt depends on {kx,ky), and 
that a gapless point exists at fc^ = at some (k^^ky) near 
the fc-point shown in the figure. In order to understand 
this behavior completely, we must perform many more 
expensive calculations. Therefore, we leave further study 
of this problem for future studies. 

Finally, we note the possible Fermi surfaces. There are 
two kinds of possible Fermi surfaces (lines) at low tem- 
peratures if there are crossing points at some (fc^, ky)^ de- 
pending on the value of kz . One possibility is Fermi lines 
stretching along the c*-axis. However, it is not natural for 
the Fermi level to always be at the crossing points. The 
other possibility is more probable: acicular Fermi sur- 
faces pointing along the c-axis. The Fermi velocity along 
the c*-axis is much smaller than those along the a6-plane. 
Note that the density of states at < 1 meV~10 K 
may be changed due to this three-dimensionality if the 
dispersion along the c*-axis is relevant. 
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Table IV. Energy difference (meV) along the c*-axls at ambient 
pressure (ap), with uniaxial strain of 2 kbar along the a-axis (a2) 
and 3 kbar along the t>-axis (b3), at 8 K and room temperature 
(RT). Values show E{k±, kz = 0) — E{k±, kz = tt) at each k± = 
{kx, ky). 



T=8 K 


width 


ap 


a2 


b3 


( — TT, 7r) 








LUMO 
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2.6 
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HOMO 
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-4.0 
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-2.0 
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3.0 
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-2.3 


-2.3 
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-0.2 


-0.2 


-0.7 


T=RT 


width 
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2.3 
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HOMO 


-3.8 


-3.7 


-3.8 
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LUMO 


-2.3 


-2.4 


-2.4 


HOMO 


2.7 


2.7 


2.8 


(-7r,0) 








LUMO 


-2.2 


-2.3 


-2.6 


HOMO 


0.2 


0.2 


0.1 



3. Summary 

We have relaxed atomic positions of a-(BEDT- 
TTF)2l3 at ambient pressure and with uniaxial strain 

along the a- and 6-axes to determine the electronic struc- 
ture by the first-principles calculation. We have con- 
firmed the existence of anisotropic Dirac cone dispersion 
in all the cases. Uniaxial strain of 3 kbar along the fo-axis 
slightly increases the band widths, while uniaxial strain 
of 2 kbar along the a-axis decreases or does not change 
the bandwidths, though the overall effects of pressure 
to the bandwidths are not large. The effects of thermal 
contraction on the electronic structure are much stronger 
than those of uniaxial strain. There are two main effects 
of temperature. One is the increase of the bandwidths 
with decreasing temperature. The other is the change 
of the electronic structure near the chemical potential. 
Only the Dirac cone dispersion appears near the chem- 
ical potential in the energy range of the temperature at 
8 K, while the electron pocket at (— 7r,7r) and then the 
hole pocket at (— tt, 0) appear within the energy range of 
the temperature, as well as the Dirac cone dispersion, at 
RT. This Dirac cone dispersion yields the linear density 



of states, which gives the dependence of the carrier 
density up to about 100 K. If the effects of uniaxial strain 
of 10 kbar along the a-axis are almost the same as those 
of 2 kbar strain, the theoretical carrier density is compa- 
rable to the experimental value. The effects of uniaxial 
strain are different between the a-axis and the 6-axis. 
uniaxial strain along the a-axis gives almost the same 
density of states as that obtained at ambient pressure, 
at least in the energy range of temperature, while strain 
along the 6-axis gives a much large value of the density 
of states at 8 K. It is apparent that the effects of uni- 
axial strain are different between the a- and b- axes on 
a low-energy scale. These theoretical findings may have 
some connections with the experimental phase diagram 
that depends on the direction of strain. 
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Fig. 2. Electronic structures at ambient pressure (ap,») and with 
uniaxial strain of 2 kbar (a2,+) along the a-axis and 3 kbar along 
the 6-axis (b3,x) at RT (a) and at 8 K (b). A comparison of RT 
(x) and 8 K (+) at ambient pressure is shown in (c). We take 
E = for the chemical potential. 




Fig. 3. (a) Electronic structure at fc^ = near anisotropic Dirac 
cone dispersion at ambient pressure and RT. Contour plot is also 
shown at the bottom. Solid (dashed) lines are for the HOMO 
(LUMO) band. \vF^y\ is much smaller for ky < kyo than for 
ky > kyo in the HOMO band, while \vF,y\ is much smaller for 
ky > kyo than for ky < kyo in the LUMO band, where kyo is 
the wavcnurnbcr along the 6*-axis at the crossing point of the 
anisotropic Dirac cone. Projections of E{kx, ky, kz = 0) in the 
area from (— tt, 0) to (0, tt) onto the plane perpendicular to the 
afe-plane at ambient pressure (ap), with uniaxial strain of 2 kbar 
along the a-axis (a2) and 3 kbar along the fc-axis (b3) (b) at room 
temperature, and (c) at 8 K . Vertical lines at (0,7r) and (— 7r,0) 
are guides for the eye. (d) Schematic representation of projection 
of E{kx,ky, 0) in (b) and (c). We do not remove hidden surfaces 
in the figures. The finest mesh is 324x324 in all cases. We take 
E = to be the crossing energy of the anisotropic Dirax; cone 
dispersion. 



J. Phys. Soc. Jpn. 



Full Paper 



Hiori Kino and Tsuyoshi Miyazaki 9 




Fig. 4. Crystal structures of conducting plane (the afe-plane) and 
transfer integrals in Kendo's crystal axes.^' ^ Molecules A and A' 
are connected by inversion, al', a3' and a4' are transfer integrals 
to the next nearest neighbor molecules along the a-axis. 
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Fig. 5. Two-dimensional band structure in the scmi-cmpirical ap- 
proximation, which is compared with Fig. 2. The scmi-cmpirical 
transfer and overlap integrals are from ref. 7. Note that we do 
not use the fitted transfer integrals from Table III in this figure. 
We show this figure because the axes and fe-points differ from 
study to study. 
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Fig. 6. The density of states at room temperature (RT) (a) and 
8 K (b) at ambient pressure (ap, solid lines), with uniaxial strain 
of 2 kbar along the a-axis (a2, thin dashed lines) and 3 kbar along 
the fe-axis (b3, thick dashed lines), (c) Effective carrier density 
estimated as n(T) = j]^^ D{E,T)dE/V{T), where D(E,T) is 
the density of states, and V(T) is the cell volume at temperature 
T. We interpolate values of the transfer integrals and the cell 
volume linearly as a function of temperature to calculate D{E, T) 
and V{T). 
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Fig. 7. Electronic structure along c*-axis near {kx,ky) of the 
crossing point of the Dirac cone dispersion. We talce the gap- 
less point to be at E = 0. 



